function y = rst_l(temp, rym)
    Fs = 8192;
    wavef = [261, 293, 329, 349, 391, 440, 494];
    x = linspace(0, 2 * pi * rym, floor(Fs * rym));
    y = sin(wavef(temp) * x) .* (1 - x / (2 * pi * rym));
end